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Abstract. The relation between the Toda lattices and similar nonlinear chains 
and orthogonal polynomials on the real line has been elaborated immensely 
for the last decades. We examine another system of differential-difference 
equations known as the Schur flow, within the framework of the theory of 
orthogonal polynomials on the unit circle. This system can be displayed in 
equivalent form as the Lax equation, and the corresponding spectral measure 
undergoes a simple transformation. The general result is illustrated on the 
modified Bessel measures on the unit circle and the long time behavior of their 
Verblunsky coefficients 



1. Introduction 



In 1975 J. Moser ^2 ^] suggested a method for solution of the finite Toda 
lattice equations (specifically, the Cauchy problem for such lattices) based on the 
spectral theory of finite Jacobi matrices. Later on Yu. M. Berezanskii [S] adapted 
this method to semi-infinite Toda lattices 
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where ' means differentiation with respect to t, in the class of bounded real b's and 
positive a's with the initial data {6 n (0) = 6„(0), a„(0) > 0}. The key idea is to 
compose a semi-infinite Jacobi matrix 
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and trace the evolution of the matrix-valued function J = J(t) and its spectral 
characteristics. It turned out that l|l.lfl can be paraphrased in equivalent form in 
terms of J itself (the Lax equation) 
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with the standard notation X± for the upper (lower) projection of a matrix X, as 
well as the corresponding spectral (orthogonality) measure dfi(x, t) which undergoes 
a simple transformation 

(1.5) dn(x,t) = e~ xt d^(x,0). 

Hence the solution of (|l.ll) boils down to a combination of the direct spectral 
problem (from {a„(0), b n (0)} to dfj,(x,0)) at t = 0, plus 1(1. 5J) . phis the inverse 
spectral problem (from d[i(x,t) to {a n (t),b n (t)}) at t > 0. 

The theory of orthogonal polynomials on the real line plays one of the first fiddles 
in the performance (albeit not entering the final result directly). For instance, it 
furnishes a nice setting for solving the inverse spectral problem. There is a parallel 
theory of orthogonal polynomials on the unit circle (OPUC) which has experienced 
a splash of activity lately thanks to primarily Simon's disquisition 120] . So the 
question arises naturally whether there exist nonlinear chains (so to say, the "Toda 
lattices for the unit circle") which can be handled by the similar method. The 
main goal of the present paper is to develop the "Moser-Berezanskii scheme for the 
unit circle" based on the spectral theory of a certain class of unitary matrices in 
application to a system of nonlinear differential-difference equations known as the 
Schur flow. 

We begin with some basics on orthogonal polynomials on the unit circle (OPUC). 
Given a nontrivial (i.e., not a finite combination of delta functions) probability 
measure /z on the unit circle T with the moments 

Hk := [ C k dfi, k£ Z= {0,±1,...}, 
Jt 

we define the monic orthogonal polynomials <& n (z, fi) (or just $„ if \i is understood) 

by 

(1.6) / C fe *«(C)^ = 0, fc = 0,...,n-l; $„(z) = z n + 1^ 1 + . . . + $„(0). 
Jt 

Clearly such system is uniquely determined and 

(1.7) / $ m (C)$^0 dn = 0, m £ n. 
Jt 

The orthonormal polynomials tp n = n n ^ n , n n > enjoy the property 

< f fm(C)<< 5 n(C) d(J< = 5 mn . 

A key role throughout the whole OPUC theory is played by the sequences of 
complex numbers {a n }n>o> \ a n\ < U 

(1.8) a n = a n (n) = -$ n+1 (0), n £ Z+, a_i := -1, 

known as the Verblunsky coefficients or parameters of OPUC system. Firstly, due 
to the celebrated Verblunsky theorem, there is one-one correspondence between 
the class V of all nontrivial probability measures on T and the set D°°, so each 
sequence of complex numbers {7 n }n>o from the open unit disk D comes up as a 
system of parameters for uniquely determined measure fi £ V . Secondly, Verblun- 
sky coefficients (11.8(1 enter the Szego recurrence relations given in the vector form 

by 
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is called the Szego matrix, and so 



%i{z) = T n (z)T n _x(z) . . . T (z) 



*n+l(«) 1 _ r M [ 1 

is the transfer matrix. So both monic orthogonal and orthonormal polynomials are 
completely determined by the sequence a n , the latter because of the equality 

n-l 

^^na-K-i 2 )- 

j=o 

To complete with the basic properties of OPUC let us mention explicit determi- 
nant formulae for both monic polynomials and Verblunsky coefficients in terms of 
the moments of the orthogonality measure: 
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One of the most 
is the discovery by 
multiplication by £ 
if |ro — ra| > A: for 
the real line case) 
by orthonormalizin 
expressed in terms 

(1-11) X 2n(z 



interesting developments in the theory of OPUC in recent years 
Cantero, Moral, and Velazquez [S] of a matrix realization for 
on i 2 (T,/i) which is of finite band size (i.e., |(CXro>Xn)| = 
some k; in this case, k — 2 to be compared with k = 1 for 
Their basis (complete, orthonormal system) {xn} is obtained 
g the sequence 1, C, C ■> C 2 j C~ 2 > ■ ■ ■• Remarkably, the x's can be 
of <^'s and ^*'s (see [03 Proposition 4.2.2]) 

) = Z~ n lf2 n (z), X2n+l(z) = Z~ n ip2n+l(z), Tl G Z+ , 



and the matrix elements 

C(M) = IKmll = (CXm,Xn), 

in terms of Verblunsky coefficients 
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V : 

< p < 1. 



There is an important relation between CMV matrices and monic orthogonal 
polynomials akin to the well-known property of orthogonal polynomials on the real 
line: 

(1.13) =det(z/„-C (n) ), 
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where C^ n ' is the principal nxn block of C (see, e.g., ^3 p. 271, formula (4.2.47)]). 

The CMV matrices C play much the same role in our study of the Schur flows 
that Jacobi matrices (|1.2|l in the Toda lattices setting. 

We are in a position now to announce our main result. 
Theorem 1. Let a n (t), n G Z + , be a sequence of complex valued functions with 
\®n(t)\ < 1 for t > and let a_i = — 1. The following three statements are 
equivalent. 

(1) a n solve the Schur flow equations 

(1.14) a' n (t) = (l-\a n (t)\ 2 )(a n+1 (t)-a n - 1 (t)), t > 0; 

(2) The CMV matrices C(t) satisfy the Lax equation 

(1.15) C'(t) = [A,C], 

where A is an upper triangular and tridiagonal matrix 



(1.16) A 



( 3?a poA papi 

— 3fJa 2 ai P2A2 

v ■■■ ■■• - j 



A n := a n+ i(t) - a n -i(t); 
(3) The orthogonality measure fi(£,t), having a n (t) as its Verblunsky coeffi- 
cients, satisfies 

(1.17) d^t) = C{t)e t ^ 1 U^Q), 

where C is a normalizing factor. 

We refer to (|1.17fl as the Bessel transformation of the measure dfi = d/i(C, 0). 

Remarks. 1. It is not hard to see that if {«„(£)} solves ()1.14f) with the initial 
data |a„(0)| < 1, then |a ra (t)| < 1 for each t > 0. Indeed, assume for the contrary 
that for some n > there is t > such that |a n (io)| = 1 an d |a„(i)| < 1 for all 
< t < t . It follows from l|1.14|) and p n = 1 - \a n \ 2 that 

(PnY = 2pnPn = - a l"n - ot n a' n = -2p 2 JR(a n +ia n - a n a n -i), 

and so 

(1.18) /4 = -p„5R(5„+ia„ - a„Q n -i), < t < t . 
Hence 

(1.19) p n (t) = p n (0)exp(- J $t(a n+1 (s)a n (s) - a n (s)a n _i(s))ds| , 

and the right hand side is bounded away from zero as t — > £0, whereas the left hand 
side tends to zero. The contradiction shows that |a n (t)| < 1 for all n E Z + and 
t > 0. Therefore, by (3) the solution of the initial-boundary-value (IBV) problem 

(1.20) a' n (t) = {l-\a n (t)\ 2 ){a n+1 (t)-a n _ 1 (t)) ) t > 0; |a„(0)| < 1, 
n € Z + , a_i = —1 exists and unique. 
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2. We can modify the second statement by observing that 



A + A* = C + C* = 
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and so 
(1.21) 



C'(t) = [B,C], 
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which makes it closer to its counerpart in the Toda lattices setting. 

So, once again, the solution of IBV problem (|1.2UI) amounts to a combination 
of the direct and inverse spectral problems with l|1.17fl m between. Note that 
the orthogonality measure /i(£,0) can be retrieved from the initial data a„(0) by 
either the Spectral Theorem for the CMV matrix C(0) (|1.12|) or via orthonormal 
polynomials, since /u(£, 0) arises as a *-weak limit of the sequence of measures 
\(p n \~ 2 d'm, dm being a normalized Lebesgue measure on T (Rakhmanov's theorem). 
In turn, the Verblunsky coefficients a n (t) are recovered from the measure n(C,t) 
by flit)! . 

The Schur flow (|1.14(l emerged in 0E] under the name discrete modified KdV 
equation, as a spatial discretization of the modified Korteweg-de Vries equation 

dtf = 6/ 2 clJ - dlf. 

In [S] the authors deal with finite real Schur flows and suggest two more distinct 
Lax equations based on the Hessenberg matrix representation of the multiplication 
operator (see also [5]). In |13U14) the Bessel modification of measures appeared and 
a part of our main result which concerns (3) (1) is proved. In a recent paper 
the author deals with the Poisson structure and Lax pairs for the Ablowits-Ladik 
systems closely related to the Schur flows. The latter can also be viewed as the 
zero-curvature equation for the Szego matrices (cf. [S]) 

r n (z, t) + T n (z, t)W n {z, t) - W n+1 (z, t)T n (z, t) = 0, 



W n (z,t) := 



+ l-a n _ia n -a n -a n -\z 
-a n -iz-a n l-a n -ia n + z~ 

We proceed as follows. The proof of our main result is presented In Sections 2 and 
3 with some comments on the general IBV problem and doubly infinite systems. We 
also show that some properties of Verblunsky coefficients for the Bessel transformed 
measures, such as the rate of decay, are inherited from those of the initial data. In 
Section 4 we study the modified Bessel measures and polynomials on the unit circle, 
a nice example which corresponds to the zero initial conditions in our setting. In 
this case the long time behavior of the Verblunsky coefficients can be obtained. 
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2. Proof of Theorem 1: (1) <^> (2). 

(1) =>■ (2). Once the Lax pair is enunciated, the proof goes through by brute force 
computation which is much more involved compared to the Toda case. 

Let 2e m := 1 - (-1)™, m € Z+, and e_i = 1, so {e m } m > = {0, 1,0, 1, . . .}, 

£m + £m + l = 1, £m £ m+l = 0, £ m — £ m +l = (— l) m+ . 

It is instructive to write the diagonals of C (|1.12f) in a unique way 

(2.1) Cmm Q; m Q; m _i, 

(2.2) c m+2,m PmPm+l^rrn c m,m+2 > 

and 
(2.3) 

(2-4) C m . m +i d m _i/3 m £ m . 

In the same vein for the matrix entries of A (|1.16(l 

(2.5) &mm ^C^mC^m— 1) &m,m+2 — PrnPm+li 

and 

(2.6) ^m,m+l — ~t~ PmA m £" m . 

Next, it follows from lfTTTl|l and (fT~T%|) that 

(pmPm+i)' = —PmPm+i^(a m+ 2a m+ i — a m a m -i), 
(a m a m _i)' = a m _ip m A m + a TO p m _ 1 A rn _i, 

and 

(a m -ip m )' = [p m _ 1 A TO _i — a m _i5R(a m+ ia m — a m a m ~i)\ p m , 
(a m+ ip m y — [p^ +1 A m+ i — a m+ i3?(a m+ ia m — S TO a TO _i)] p TO . 
Hence, for derivatives of the CMV matrix entries we have now 

(2.7) c' mm = -(5 m a m -i)' = -a m -i/»„A m - amp^jAm-i, 

(2.8) c ' m ,m+2 = — PmPm+l$l(®m+2(Xm+l ~ ~&rn a m-l)£m+li 

(2.9) c' mm _ 2 = —Pm-2Pm-l$l((Xm tm-l — ®m-2(Xrn-a)£rn, 
and 

C 'm,m+1 = -ara+lS(a m 4-ia m -5 m « m _i)] p m £ m +l 

(2.10) - [p m _ 1 A m _ 1 - a m _i^(a m+ ia m - a m a m _i)] p m £ m , 

c m,m-l = [^Am -5 m S(a m a m -l -Q m -lQ ra -2)] p m -l£m+l 

(2.11) - [ / 9^ l _ 2 A„ l -2 — a m _ 2 ^(a 

Let if = H-ffronll = [A, C], if mn = X)i( a ™ c ™ ~ c mi a in ). Since both A and C are 
of band size 2, i.e., a mn = c mn = for |?ti — n| > 2, and A is an upper triangular, 
i.e., a mn = for m > n, we actually have 

(2.12) K mn drnmCmn ~l~ ^m,m+lCm+l,n "I - ^m,m+2Cm+2,n 

CmnG'nn Cm,n — l^n — l,n ^^-2^-2^' 

We want to show that 

(2.13) (t)=K j(t), m,m + j£ Z + . 
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To this end we will plug (2.1)-(2.6) into (|2.12(1 and compare the outcome with 
(2.7)-(2.11). 1 For j > 5 and j < —3 the equality holds for trivial reason, as the 
both sides in (|2.13|) vanish. For j = 4 

-^m,m+4 ^m,m+2Cm+2,m+4 ^m,m+2^m+2,m+4 

= pmPm+lPm+2Pm+3{£m+S ~ £m+l) = 0. 

For j = 3 by (2.2) and (2.4)-(2.6) 

^m,m+3 ^m,m+lCm+l,m+3 ~l~ ^m,m+2^m+2,m+3 

PmPm+lPm+2 [£m(A m — Ck m +i + a m _i) + £ m +l(o ; m+3 — a m+l — A m+2 )] = 0, 

that is consistent with the banded structure of size 2 of C 
The main work begins when \j\ = 0, 1, 2. 

1. j = -2. We have by (2.2) and (2.5) 

(2-14) K mrn — 2 — £m t 7n— 2{p , rnm ^m— 2,m— 2) 

= pm-2pm-ffi(&m-2Q!m-3 ~ S m a m _i)e m , 

and so (f^7T5f> holds by (2.9). 

2. j = -1. Write -Km,m-i = xgL-i + ^ 2) m -i with 

-^m m— 1 Cm,m—l( y Q , rnm ^m— l.m— l) 

= (a m Pm~l£m+l — Q!m-2Pm-lfm)!ft(ain-lQ!m-2 — a m a m _l), 

mm— 1 ^m,m+l^m+l.m— 1 ^m— 2,m — lCm,m — 2 

= Pm—l(P m Am£m+l ~ ftn-2^m-2 £ m) 

and hence 

(2.15) if mjm _i = [p m A m — a m $t(a m a m -i — a m _ 1 Q! m _2)] p m _i£ m+ i 

_ [pL-2^m-2 - «m-2S(«m«m-l ~ a m _J.a m _2)] p m -lE m , 

Now (|2~T3ll follows from (2.11). 

3. j = 0. Write 2f mm = Kmm + ^mm with 

-^mrn — ^m,m+2^m4-2,m Cfn,m- 2^"m— 2,m — [(PmPm+l) (Pm— 2pm— l) ] ^m? 
(2) 

mm — ^m,m+lCm+l,m Crn,m— l^m— l,m 

= (S m -|-lp m A m + Q m _2ft 2 „_iA m _i)£ ra 

^ (ckm-lPm^m + a mP m -l ^m- 1 )£m+l ■ 

But 

(PraPm+l) 2 ^ (Pm-2/5m-l) 2 + Ot m +ip m A m + a m _2Pm-l ^m-l = 
Pm(l ~"m+l«m-l) ~ P m -l(l -amOm-2) = -"m-lpiAm — O-mPm- 1 ^m- 1 , 

and so by (2.7) 

K-mm (f^m— iPjn^'m ^mPm-l^™-l)(^™ £"m+l) — ^mm' 



^One needs some paper and patience to slug one's way through the lengthy calculation. 
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4. j = 1. Now K m , m+1 = K^ m+1 + K%] m+1 + K^ m+1 with 

m+1 ^m,m+l(^mm ^m+l,m+l) 

— (^m+l^m+1 — ^m-l^m)/'m^(^m+l^m G^m^m — l); 
m,m+l ^m,m+l l,Cm+l,m+l ^rrim ) 

(^m^m+1 ^m^m)Pm(^m+l^m ^m^m — l); 
m m+1 ^"m,rn-\-2^m-\-2.m-\-l ®"m — 1 , m + 1 ^-m , m — 1 

Hence / o~ 1 i^ m , m+ i = u m e m . + w m £ m+ i with 

u m = -a m -iK(a m +ia m — a m a m _i) + ^m-2p m _i — a m p^ +1 



^m (^m+l^m ^m^m — 1 ) j 

) + r m— 1 

A m (a m+ ia m - a m a m _i). 



Next. 



u (l) := «_ - -^_„ 2 



m—l ^m+l) (^m+l^m C^m^m — l) 

a m +ia m a m _i — Smtt 2 ,^! +a m+ iQ! TO Q; m _i — a m + a m _2P rn _ 1 

m — l ^m 

m — l ' 



and so 



In exactly the same way 

v m = -a m+1 dl(a m+1 a m - a m a m -i) + p m+ tA m+ t, 

and finally 

(2.16) K„^ m+ i = [p 2 ra+1 A m+ i - a m+1 dl(a m+1 a m - a m a m -i)] p m e m +i 

~ [Pm-i^m-i - %-iK(5ro+iai m - a m a m _i)] p m £m- 
We come to (I2.13|) on account of (2.10). 

5. j = 2. Write tf m , m+2 = K$ m+2 + K% m+2 + K^ m+2 with 

-^m,m+2 = a m,m+2{c m +2,m+2 ~ C mm ) — Pmpm+l( a m a m-l ~ &m+2&m+l){£m + £m 

-^m,m+2 = c m,m+2( a mm ~ a m+2,m+2) = PmPm+l£m+l^(~(* m a m -i + a m+ 2<^m+l), 
7^(3) _ 

m m+2 — ^m,m+lCm+l,m+2 Cm,m+l a m+l,m+2 

= PmPm+i [(«m+2A m + a m -iA ra+ i)e ra — (a m A m + a m +i A m+ i)e m+ i] . 

We have 

(* m + 2 A m + a m -iA m+ i = — a TO a m _i + a m +2a m +i 
draAm + a m+ iA m+ i = — a m a m -i + a m+ 20i m +i 
and it follows from (2.8) that 

(2-17) K m m+ 2 = — PmPm+l^(ctm+2Ctm+l — CLm&m- 1 )£m+l , 

so (12.1311 holds again. The proof is complete. 
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(2) => (1). The problem we are faced with is that, in contrast to the Toda lattices, 
no a's in a pure form appear among the matrix entries of C. 
Write \a n+ ip n \ 2 = (1 - p 2 n+1 )p 2 n =(% - (p n p n +i) 2 and so 

(PnY = 2pnPn+l(pnPn+l)' + {a n +lPn)' (O-n+lPn) + (p-n+lpn) ("n+lPn)'- 

Hence 
(2.18) 

Pn — Pn+1 (pnPn+i)' + ~3l[a n+1 (a n+1 p n y]. 
Next, (a n -ip n )' = a' n ^ x p n + a n -ip' n , so that 



(2.19) 



<*n-l = [(an-lPn)' ~ On-lPn] ■ 

Pn 



The right hand side of (|2~T3|) can be expressed in terms of derivatives of the 
CMV matrix entries and thereby, via the Lax equation, of a's themselves. First, 
by ifTTIjl and (|?T7|) 

(2-20) (pnPn+l)' = dn,n+2 + c 'n+2,n = K n,n+2 + K n +2,n 

(2.21) = p n p n+ idi(a n an-i ~a n +20!n+i)- 

Next, it is immediate from (2.3), (2.4) that a n +\p n = c n+ i tn e n + c n! „+i£„+i, and 
so by (215) . (|2T^|) 



c n+l,n e n + c n,n+l £ n+l 



K r 



K 



n,n+l^"n+l 



cXn+iPn^(a n a n -i - a„+ia n ) + p n p 2 l+1 A n+1 



(2.22) (a n+1 p n )' 
(2.23) 
Similarly 

{otn-xpn)' = oi n -ip n U(a n a n -i — a n+ ia n ) + p n p n _ x A n _i. 
Plugging fFZfy and $TF^i into gives 

Pn = PnP n +i$l(®nOin-i - a n+2 a n+1 + a„+iA n+ i) + p n \a n+1 | 2 3f?(a„a„_i - a n+1 a n ) 
= PnP n+ i^{-a n+ 2an+i + a n+ ia n + a„+i A n+ i) + p„3?(a„a„_i - a n+ ia n ) 
= p n 3?(a n a„_i - a n+ ia n ). 

In the upshot, the Schur flow equations emerge from (|2.19(l : 

ol n _xPn = a„_ip n 3?(a„a„-i - Sn+icnn) + pnP„-iA n _i 

= -Ot n -lPn$l(Un(Xn-l ~ Un+l&n) = PnP n -l&n-l; 

as claimed. □ 
Remark. We could equally well have considered the general IBV problem, that, 
strictly speaking has nothing to do with OPUC: 

o'Jt) = (1 - \a n {t)\ 2 ){a n+l {t) - a„_i(t)), t > 

with a continuous boundary function < 1. The above evaluation shows 

that the Lax form of such problem is C' g = [A g ,C g ] with 

\ 



c 9 = 





&iPo 


A) Pi 








a-iPo 


-aiao 


-a pi 











a 2 pi 


— a 2 ai 




P2P3 







—a\p2 


— 0.30.2 


-a 2 p3 











O4P3 


-0^03 
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and 

( -Kaoa-i poAo popi 

— 5Ra 2 ai p 2 A 2 P2P3 



T„ 



v '•• ••• / 

(to be compared with (|1.12() and (|1.16(l ). Furthermore, the doubly infinite system 

a' n (t) = (1 - \a n (t)\ 2 )(a n+1 (t) - a„_i(i)), t > 0, neZ 
with the initial data {a„(0)}fe e z C D has its equivalent Lax form 

C' = [i,C], 

where C and A are doubly infinite extensions of l|1.12|l and l|1.16(l given by the same 
expressions (2.1)-(2.4) and (2.5)-(2.6), respectively, with 

1 - (-l)H 
£ m = ^ ' m€Z. 

3. Proof of Theorem 1: (2) => (3) (1). 

(2) => (3). Let 

fl,(t) := (C - ziy 1 = \\r mn (z,t)\\™ n=Q 

be a resolvent of the CMV matrix (|1.12|l . It is easy to see that R z obeys the 
same Lax equation l|1.15fl . Indeed, differentiating the identity (C — zI)R z — I with 
respect to t entails 

C'{t)R z {t) + (C(t) - zI)R' z {t) = 0, 

and so 

(3.1) R' z (t) = -R z (t)C'(t)R z (t) = -R z (t)[A,C]R z (t) = [A,R Z ], 

as claimed. 

Take the equation for (0, 0)-entry of l|3.1|l : 

r' 00 = $ta Q r QQ + p A Q ri Q + p pir 2 o - r o5iao 

= p (ai + l)n + PqP\T20- 

As it follows from (C — zI)R z — I 

(So ~ z)r Q0 + aiPo^io + PoPir20 = 1, 

which allows to eliminate r 2 o in favor of rio, ^oo, and so 

(3.2) r' Q0 = p Q r 1Q + 1 - (a - z)r 00 . 

By the Spectral Theorem, the resolvent entries can be found from 



r m , n {z,t)= / dfj,{(,t) 

with x„ which is particularly simple for the first two elements 



rooOM) = / — : , r 10 [z,t)= / dp,(C,t), 



^i(z) = p 1 (z — ao). Hence for the right hand side of (|3.2() 

C-2 



, C + C - 23ffa 
Ajrio + 1 - (ao - 2)^00 = / t ap((,t) 
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holds, and we end up with a differential equation for orthogonality measures 
V(C,i) = (C + C-2Rao)d/i(C,t). 

Finally, 



d/i(C,t) = exp|y (C + C-23?a (s))ds| <MC,0) 

as needed. 

(3) => (1). We start out from the transformation of the orthogonality measure 
d/i(C, i) = <p(C; d/z(C, 0) and derive a differential equation for the moments 



Now 



and so 



exp(^(C + C" 1 )) 
/(*) 



/(*) = / e 



,t(c+c _1 ) 



dM(C,0), 



/'(*). 



y'(C, *) = (C + r X MC, t) - 7^yv(C *)• 

As J'/ -1 = ci + c_i, we have 

(3.3) /4(t) = Mfc-i (*) + Mfc+iW - 9(t)nk(t), g(t) = ci(t) + c_i(t), fc e Z. 

The rest is based heavily on (|1.1(J|I which relates Verblunsky coefficients and 
moments of the orthogonality measure. The idea to differentiate determinants and 
take into account l|3.3|) goes back to |S] , see also ^5] Lemma 1] . For a set of integers 
k\ < &2 < ■ ■ ■ < k n denote 



T{k\ , . . . , k n 



( Mfci Mfei-i 
fj-k 2 Mfc 2 -i 



Mfcl-n+1 \ 
Mfc 2 -n+l 



\ Mfc„ Mfc„-1 •■■ Mfc„~n+1 / 



D(kx,...,k n ) := detT(fci,...,fe n ), and so L>„ = D(0, 1, . . . , n - 1). Put G r , 
D(-l, 0, . . . , n - 2) and write fl~TU|) as $„(0, <) = (-1)"G„L>- 1 . Then 



(3.4) 



r" n —CD' 



It is clear from (|3.3(l that intermediate determinants in the sum 



DL = 



Mo M-i 
Mi Mo 

Mn-1 Mn-2 



M-n+1 
M-n+2 

Mo 



Mo M-i 
Mi Mo 



Mn-1 Mn-2 



M-n+1 
M-n+2 



Mo 
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have the same value —gD n , whereas the first and the last ones equal, respectively, 



Mo 
Mi 



m'-i 
Mo 



Mn-l Mn-2 



Mo 
Mi 



M-i 
Mo 



M-n+1 
M-n+2 



Mo 

M-n+1 
M-n+2 



Mn-l Mn-2 



Mo 



£>(-!, 1, 2,..., n -!)-<?(*)£>„, 



= D(0,l,...,n-2,n)- fl (t)D n . 



Hence 



(3.5) = -ng(t)D n (t) + D(-l, 1, . . . , n - 2, n. - 1) + L>(0, 1, . . . , n - 2, ra). 
Similarly, 

(3.6) G^(i) = -ng(t)G n (t) + D(-2,0, . . . ,n-3,n-2) + D(-l, 0, . . .,n-3,n-l). 
After plugging Ij3.5|l and 13.6|l into (|3.4(1 we come to 



f-lV 

<(0,<) = ^-^{^[^(-2,0, 



,n - 2) + £)(-!, ■ • • ,n - 3,n - 1) 



- G n [£>(-l,l,...,n-l) + I>(0,...,n-2,n)]}. 
Let us now go over to the right hand side of Q1.14J1 . written for & n (0): 



(i - |$„(o)| 2 ) ($ n+1 (o) - $^(0)) = (-1)"+ ^ f"' 2 {^±i 



G n -i 
D n -i 



The standard Silvester identity applied to the matrix T(0, 1, . . . , n) gives 

D„ — \G n \ 2 = D n+ iD n -i, 

and so 

(-l) n+1 

(1 - |$„(0)| 2 ) ($„+i(0) - $„_i(0)) = 1 ' (C? n+ i£> w -i - A i+ iG„-i). 

Another application of the Silvester identity (in a bit modified form) shows that 
G n+1 D n -x = G n D(-l,l,...,n-l)-D n D(-2,0,...,n-2) 
£>„+iG n _ a = D n D{-l,0,...,n-3,n-l)-G n D(0,l,...,n-2,n), 

and we arrive at the Schur flow (|1.14|) . That completes the proof of Theorem 1. □ 
There is yet another way to prove (3) =>• (1), which gives not only l|1.14f> . but 
the differential equations for the monic orthogonal polynomials. I learned it from 
QUI Section 8.3]. 

Theorem 2. The monic polynomials $„(-,i) orthogonal with respect to the Bessel 
transformation ( |j.J7| | satisfy the differential equation 

(3.7) <(M) = $n+i(z,t) - (z + a n a n -i)^„(z,t) - (1 - |a„-i| 2 )$„-i(z, f). 

Proof. The idea is to differentiate orthogonality relations l|1.7|l with respect to t. 
Now «/ = (£ + Qw, an d we have for m < n 



(3.8) 
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It is clear that $' m is a polynomial of degree at most m — 1, and so the first integral 
in the above sum is zero. If m < n — 2 then 



f * m (*{, + C*„) <*/*(£*) =0, 
JT 



and hence 

(3.9) $;o) + *$ n ( 

with some parameters x n , y n , z n depending on t. By matching the coefficients for 
z n+1 and z n in (|3.9|) and using l|1.6(l we find 

X n = 1) 2/n = in(*) - Wl(*)- 

To get z„ take (13.81) with m = n — 1 and apply the Szego recurrences 9fl : 

0= / $„_i(^ + c$„)rfMc,*)+ / 4Kc,*), 

JT JT 

f + C#„)d/i(C,*) = -||*n||J, 

JT 

4-1 | ,2 -i 

— 5— = On-1 - 1- 



and so 



~ ||$ < J 
To find the expression for l n we turn to (|1.13|) 



n—i 

$ n (z) = z n + InZ* 1 - 1 + ...= Y[(z + a a - x ) + ..., 

3=0 

so that l n = Y^=o <*j a j-ii an d we come to l|3.7|) . Putting z = yields (|1.14[l . as 
was to be proved. □ 

It might be worth pointing out that some properties of Verblunsky coefficients 
for the Bessel transformed measures (such as the rate of decay) are inherited from 
those of the initial data. 

Theorem 3. Let a n (0) enjoy either of the properties 

(1) {a n (0)}e£P, p=l,2; 

(2) |a„(0)| < Ce- an , a>0. 

Then the same holds for a n (t) for each t > 0. 

Proof. It is obvious from (| 1 . 1 Tf> that [i(C,t) belongs to the Szego class (i.e., log^' £ 
A 1 (T) if and only if 0) does, and so the first statement with p = 2 follows from 
fundamental Szego's Theorem |191 Theorem 2.3.1]. As for the case p = l, note that 
by Baxter's theorem (see, e.g., JHI Theorem 5.2.1]) //(£, 0) = wdm with w > and 
w S W, class of absolutely convergent Fourier series. It is clear from l|1.17[) and the 
Wiener-Levy theorem that 

d/x(C, t) = te(C, t)dm, w{(, t) = e t( - c+r ^w(C) e W, 

ty(C, i) > 0, and so the repeated application of Baxter's theorem does the job. 
To prove the second statement, we introduce the Szego function 



log •&?(£) <im I , 
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defined for an arbitrary measure /i = wdm + fj, s from the Szego class. A straight- 
forward computation gives for the Bessel transformed measures (I1.17|) 

(3.10) ^(z, M (i)) = ex P Q^^[i(C + r 1 )+log^(C)]rfmj = e tz D(z, f,(0)). 

Denote by Rt the radius of convergence of the Taylor series for D^ 1 (z, /J,(t)) about 
the origin. By the Nevai-Totik theorem [*H)] 

limsup |a ra (i)| ™ — R^ 1 , 

n — >oo 

and it is clear from H3.1()|) that Rt — Rq for alH > 0, as claimed. □ 
Note that under assumptions of Theorem 3 the series 

oo 

JC = -^2aj(t)aj-i(t) 

converges absolutely and the Schur flow can be written in the form 

aj. = {a j} H}, H = -2QK, 
where the (formal) Poisson brackets are defined by 

df dg df dg 



doij daj daj daj 



So l|1.14|l is the evolution of the Verblunsky coefficients under the flow generated 
by the Hamiltonian — (cf. [15|). 

4. Modified Bessel measures on the unit circle 

Because of the boundary condition a_i = —1 IBV problem (|1.20(1 with zero 
initial conditions 

a (0) = ai(0) = . .. = 
has a nontrivial solution. Theorem 1, (3), shows that we are dealing now with the 
Bessel transformation of the Lebesgue measure 

dp(C,t) = C(t)e t(c+rl) dm, 

called in the sequel the modified Bessel measures on the unit circle, with (3 n {t) the 
Verblunsky coefficients of fj,(-,t). The corresponding system of orthogonal polyno- 
mials has arisen from studies of the length of longest increasing subsequences of 
random words 0] and matrix models JSj ( see 1X01 example 8.3.4] for more detail 
about the modified Bessel OPUC). 

Note first that C{t) can be easily computed 

c -'« = -k C ht^ Cw 

71 = 







00 t 
7Z 



In 

Wt), 



^ (n!) 2 

n=0 v ; 

where Ik is the modified Bessel function of order k. Similarly, for the moments of 
the measure we have 

f I (2*) 

Jt 1 o( zt ) 
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There is an important feature of the modified Bessel measures, namely, their 
Vcrblunsky coefficients satisfy a nonlinear recurrence relation known as the discrete 
Painleve II equation (see |1(J1 lemma 8.3.5]) 2 

(4-1) 1) ^ = f3 n+1 (t) +f3 n -t(t), neZ+, 

with 0_i = -1, [3 Q = I 1 (2t)/I Q (2t). Clearly, all /3's are real now. 

Before studying the long time behavior of (3 n we prove an auxiliary result con- 
cerning the general Schur flows. 

Lemma 4. For the solution a n of the Schur flow \1. 14\ the limit relations hold 
(4.2) lim Pn (t)p n+1 (t) = 0, lim Pn (t)A n (t) = 0, n e Z+ 

t — too t — >oo 

and 



(4.3) lim lka n {t)a n ^i(t) = 7„, n e Z + . 

t — >oc 

Furthermore, "f n is monotonically increasing: — 1<7o<7i<---<1- 

Proof. Let us focus on the form of the Lax equation given in l|1.21f) - H1.22|) . As 

B* = —B, it is clear that 

(4.4) L' = {7t(L),L], L:=^-, rr(L) = L+ - L_. 

The latter is exactly what is called in [7] the Toda flow. By Proposition 5] L{t) 
converges strongly to a diagonal operator diag (do, d\, . . .), which implies (|4.2|l - 
(14. 3|) in view of the explicit expression for L. Note that (|4.3|l with j n — —d n comes 
from the diagonal entries of L, whereas both relations in 14.2fl from the off diagonal 
entries. 

The second statement follows from (|1.19|) 

Pn(t) = /5 rl (0)exp |y $t(a n (s)a n -i(s) - a n+1 (s)a n (s))ds 

and since the left-hand side is bounded as t — > oo we have 7„ < J n +i- ^ 
Let us go back to the modified Bessel measures and their Verblunsky coefficients, 
and prove 

Theorem 5. The limit relations 

(4.5) lim I3 n {t) = (-1)", 

t—*oo 



(4.6) limf(l-^(t)) = i^ 

hold for all n G Z + . 

Proo/. Note first that (gSJ is an immediate consequence of (|4.5|l and (|4.1(l . To 
prove (|4.5(l we proceed in two steps. 

1. As we know (see Lemma 4) the sequence 7„ | . Let us show that in fact all 7 n 's 
are the same. Assume on the contrary that < 7fe+i for some k. It follows from 
(p~T^i that 



2 In the notation of r n (t) = -/3„_i(t/2). 
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But then 



\Pn(t)\ V^l) » TRW St 



i(l -$(«)) tpl(t) ~ Ct 

that contradicts (|4.1|l . since the right-hand side there is bounded. So we need only 
find the common value of 7 n . 

2. We show that 7„ = — 1 by computing lim^oo Po(t). As is well known, 
I k (t) = -?L(l + o(~)), i^oo, 



J, 
and so 

lim P (t) = lim = i, 

that is, 70 = — 1, as claimed. The desired result follows from (|4.3|l by induction. □ 
Remark. It might be a challenging problem to give a direct proof of l|4.5[) based on 
the explicit formula (|1.1U|) . which now takes on the form 

Mi) = ( i)» det ll /fc ^~^ 2t ^ ^ fc '^" , neZ+, 
det ||/fe_j(2t)||o<fe,j<Ti 

and the complete asymptotic series expansion for the modified Bessel function (see, 
e.g., [13 Chapter 7.8]) 

* fv ir - (4fc 2 -l 2 )...(4fc 2 -(2 J -l) 2 ) i t 




I managed to carry out the computation for n = 1, and it seems like one needs 
n + 1 terms of this series for /?„ . 
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